Superconductivity in the repulsive Hubbard model: an asymptotically exact 

weak-coupling solution 
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We study the phase diagram of the Hubbard model in the limit where U, the onsite repulsive 
interaction, is much smaller than the bandwidth. We present an asymptotically exact expression for 
T c , the superconducting transition temperature, in terms of the correlation functions of the non- 
interacting system which is valid for arbitrary densities so long as the interactions are sufficiently 
small. Our strategy for computing T c involves first integrating out all degrees of freedom having 
energy higher than an unphysical initial cutoff f2o. Then, the renormalization group (RG) flows 
of the resulting effective action are computed and T c is obtained by determining the scale below 
which the RG flows in the Cooper channel diverge. We prove that T c is independent of flo- Using 
this method, we find a variety of unconventional superconducting ground states in two and three 
dimensional lattice systems and present explicit results for T c and pairing symmetries as a function 
of the electron concentration. 

PACS numbers: 74.20.-z, 74.20.Mn, 74.20. Rp, 74.72.-h 



I. INTRODUCTION 

The Hubbard model is widely studied as the paradig- 
matic model of strongly correlated electrons^. However, 
in more than one dimension (ID) there is controversy 
concerning even the basics of the phase diagram of the 
model. Most theoretical work on the model has focused 
on intermediate to strong interactions, U ~ W, since this 
is the physically relevant range of parameters for any of 
the intended applications of the model to real solid state 
systems. (Here, U is the repulsion between two electrons 
on the same site, and W is the bandwidth in the limit 
U = 0.) However, for such strong interactions the only 
well controlled solutions are numerical and the applica- 
tion of determinental quantum Monte Carlo methods 3 
and the Dcnsity-Matrix-Rcnormalization-Group^ have 
been limited by the fermion sigr£ and two-dimensional 
entanglement problems^ respectively. 

Here, we study the limit of weak interactions, U/t — > 0, 
where we compute the phase diagram and obtain ex- 
pressions for the critical temperatures which, assuming 
the validity of certain assumptions discussed below, are 
asymptotically exact. To be explicit, we consider the 
Hubbard model 
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for a variety of lattice systems in two and three dimen- 
sions. Here, c\ a creates an electron with spin polariza- 
tion a on lattice site i, and < i,j > and signify, 
respectively, pairs of nearest-neighbor and next-nearest- 
neighbor sites. 

Since the Cooper instability is the only generic insta- 
bility of a Fermi liquid, except for certain fine tuned 



values of t' /t and the electron density n, the only or- 
dered states that can be stabilized by weak interactions 
are superconducting states. For repulsive interactions, 
W S> U > 0, the superconducting transition tempera- 
ture has an asymptotic expansion 

T C ~W exp {-a 2 (t/U) 2 - cti(t/U) - a } 



x \ + 0{U/t) 
Wexp{-l/\pV eS ]}[l + 0(U/t)] 
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where a„ are dimensionless functions of t'/i, n and p is 
the density of states at the Fermi energy. The principal 
result we report here is to give an explicit prescription 
for computing a 2 and ol\ as a function of the electron 
density, n, and the "band structure" . On the basis of the 
present analysis, we conclude that the resulting phase 
diagram is asymptotically exact in the sense that 
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We will also explain why we are unable to give a pre- 
scription for computing o,q. In the process of computing 
ct2, one determines the symmetry of the superconducting 
ground state (e.g. s-wave, p-wave, d-wave, etc.) and the 
form of the pair wavefunction. 

There arc, of course, special situations in which a vari- 
ety of different non-superconducting ordered phases oc- 
cur. While these situations are potentially significant in 
what they imply about the behavior of the system at in- 
termediate U, in the small U limit they always involve 
a large degree of fine tuning of parameters. The canoni- 
cal example is the case of a square lattice, in which the 
model with t' = has a non-generic particle-hole sym- 
metry which leads to perfect nesting of the Fermi surface 
when the mean electron density per site is n = 1, where 
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These special situations arc thus, in some sense, not re- 
ally a part of the weak coupling problem, but rather a 
piece of the strong correlation problem that persists to 
weak coupling. If one's principal interest^— is in extrap- 
olating well controlled weak-coupling calculations to the 
range of strong interactions, this degree of fine tuning of 
the bandstructure is a small price to pay to gain access 
to phases which have broad ranges of stability for inter- 
mediate to large U . However, if we focus on the small U 
limit in its own right, then for any fixed, non-zero value 
of if /t no antiferromagnetic insulating phase occurs. 

An interesting interplay between various possible or- 
dered phases can also occur when the Fermi energy is 
coincident with a van-Hove singularity^—. While these 
singularities occur for generic values of t'/t, the density 
must be fine tuned in order for the singularity to lie suf- 
ficiently close to the Fermi energy to matter. The study 
of the behavior of the system in weak coupling tuned 
near a vanHove singularity has been explored by several 
authors^— , again as a route to understanding strong- 
coupling phases and the interplay between phases in a 
regime of parameters in which perturbativc rcnormal- 
ization group (RG) methods can be applied. However, 
again, as we are focussing on the physics of a system 
with small U /t, we do not treat the interplay with non- 
superconducting orders. 

The original idea of obtaining superconductivity from 
repulsive interactions dates back to the pioneering work 
of Kohn and Luttingcr— who derived an effective at- 
tractive interaction from the Friedel oscillations of a 3D 
electron gas. RPA calculations for a repulsive Hub- 
bard model found that near a spin-density-wave insta- 
bility there was an effective interaction which favored d- 
wave superconductivity^. The treatment of this problem 
from the standpoint of the renormalization group was 
presented by Zanchi and Schula^^ and othersii^— . 
Furthermore, the problem of competing instabilities 
of electronic systems has extensively been studied via 
the numerical functional renormalization group (FRG) 
method a 18 i 28 ~ — . While these works have made signifi- 
cant progress in our understanding of superconduvtivity 
from repulsive interactions, we present here an asymp- 
totic analysis of the problem and show explicitly the way 
the final expressions for the superconducting transition 
temperature arc independent of the initial choice of cut- 
off. Our analysis is based on the renormalization group 
framework established by Shankar— and Polchinsky22. 
In section IVIIi we present a more complete discussion 
of the relationship between the work presented here and 
previous analyses of this problem. 

This paper is organized as follows. In section II, we 
present the results of our perturbative RG treatment. 
Expressions for the quantity ct2 are presented in the spin 
singlet and triplet pairing channels. These results are 
then applied to a variety of systems in section III. Both 
lattice and continuum systems in d — 2 and d = 3 arc 
considered. In section IV, wc present the overall strategy 
of our RG calculations. In section V, wc describe the 



technical aspects of the perturbative renormalization of 
the effective interactions in the Cooper channel, and dis- 
cuss the one-loop RG flows of the effective pairing vertex 
in section VI. In section |VH] we discuss the relation of the 
present work to previous closely related approaches, and 
we conclude, in section IVIIII with a few remarks about 
future directions. There is an appendix with technical 
details. 

To simplify our notation, we henceforth adopt units in 
which h = 1, and the volume of the unit cell v = 1. 



II. RESULTS 

Before discussing the derivation, we articulate the final 
results of our analysis, i.e. we give an explicit method to 
compute cti exactly. (It is more complicated to compute 
ai, and since it is subdominant; we defer that part of the 
discussion to the more technical portion of the paper.) 
The results are expressed in terms of x(fc), the static 
susceptibility of the non-interacting system, 



X(k) 
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evaluated in the limit T — > 0, and 
d d q 



(k), 
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the density of states at the Fermi energy Ef. Here is 
the band-dispersion measured relative to the Fermi en- 
ergy, Ep, and d = 2 or 3 denotes the number of spatial 
dimensions. The integrals run over the appropriate first 
Brillouin. In general, x must be computed numerically, 
but this is a straightforward computation. 

For simplicity, we shall restrict our attention in this 
paper to systems which possess inversion symmetry and 
no spin-orbit coupling. This enables us to classify the 
possible superconducting states as having either even or 
odd parity, the former class consisting of spin singlet and 
the latter of spin triplet states. Among the even parity 
spin-singlet superconducting states, the s-wave states are 
those which transform trivially under a point group op- 
eration of the crystal, whereas the "d-wave" and higher 
angular momentum channel gap functions transform ac- 
cording to a non-trivial representation of the point group. 
The spin triplet states include p-wave and f-wave gap 
functions, all of which transform according to a non- 
trivial representation of the point group. Depending 
on the crystalline point group, the superconducting gap 
functions may transform as either a non-degenerate or an 
n-fold degenerate irreducible representation of the point 
group operations. We will present expressions for ai for 
each of these states; for any given bandstructure and elec- 
tron concentration, the physical low temperature phase 
is that one which produces the smallest value of «2- 

In order to compute 012, we must solve, for each sym- 
metry class, the eigenvalue problem displayed below. In 



2 



the spin-singlct channel, the eigenvalue problem corre- 
sponds to the integral equation 



g y k,q ^s,q n ^ s ,k 
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v F (k) 



X(k + q) +c x 



Vp 

v F (q) 
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where c\ is of order t/U, k designates a vector on the 
unperturbed Fermi surface, Sp = J dp is the "area" of 
the Fermi surface, v F (p) is the magnitude of the Fermi 
velocity at position p, and the norm of the Fermi velocity 
is defined according to 



1 

v F 



dp_ f 1 

S F \v F (p) 
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The only effect of c\ is to penalize the trivial s-wave state 
due to the bare onsite repulsive interaction. Eigenfunc- 
tions with higher angular momentum, such as the d-wave 
states, and appropriate "extended s-wavc" states are un- 
affected by ifc22r^. In the spin-triplet channel, the eigen- 
states obey 



dq 
S7 



g k,q V t,q 



1 k,q 
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Assuming that there is at least one negative eigenvalue 
present in cither Eq. [7] or [51 the quantity a 2 is obtained 
by the relation 



a<x 
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where Ao is the most negative eigenvalue. The zero- 
temperature gap function is proportional to if): 



A s(t) (fc)~T c 



v f 



^s(t)(k). 



(11) 



The computation of u\ involves an analysis of less singu- 
lar scattering processes in perturbation theory and will 
be addressed in section VI. 



III. APPLICATION TO VARIOUS SYSTEMS 

The expressions in Eqs. [Jj and [9] hold for the Hubbard 
model and are independent of the microscopic details of 
the electronic structure. In this section we determine the 
superconducting ground states for a variety of inversion 
and spin-rotationally symmetric systems in d = 2 and 
d= 3. 

There are two distinct effects of the band-structure 
which affect the asymptotic behavior of T c : Firstly, the 
existence of a superconducting instability from repulsive 



interactions at all derives from the fc-space structure of 
the effective interactions, and these depend strongly on 
n and the details of the band structure. Secondly, the 
dimensionless density of states at the Fermi energy, pt, 
varies (in 2d especially) with distance from a van-Hove 
point. To distinguish these two effects, in the figures we 
re-express the leading-order asymptotics as 



where 



T c ~exp{-[pV/ cff ]- 1 } 



V eS (n) = \X\/p = (U/t) 3 [a 2 p\ 



(12) 
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A. Rotationally invariant systems in d = 2 and 3 

We begin by considering electrons in the continuum 
limit with quadratic dispersion (which is achieved in a 
lattice system in the limit n«l): 



In 3 dimensions, 

x(q) = P 
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q - 2kf 
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where p = mkf/2-K 2 . For such a rotationally-invariant 
system, each eigenfunction can be classified according to 
its angular momentum i, and the solutions are (2£ + 1)- 
fold degenerate. The expansion of x(<z) m powers of 
q/2kf is justified when seeking the lowest angular mo- 
mentum pairing solutions and it shall suffice to com- 
pare the s-wave and p-wave solutions which require only 
the leading order expansion about q = 0. The non- 
degenerate s-wave eigenfunction has a uniform gap ev- 
erywhere on the Fermi surface. From Eq. [7J it can easily 
be seen that since the susceptibility is a positive-definite 
quantity, a negative eigenvalue for an s-wave gap function 
does not exist. However, the 3-fold degenerate p-wave 
channel with e.g. ipt(6) oc cos#, does have a negative 
eigenvalue, as can be seen by applying Eq. [9] 



MM) = \MV) 



The solution is A = — p/16 < 0. Thus, for a system with 
a quadratic dispersion and a spherical Fermi surface, the 
ground state is a p-wave superconductor to order U 2 . 
Physically, this can be understood from the fact that the 
dominant spin fluctuations are ferromagnetic. Indeed, 
such considerations have been applied extensively to su- 
perfluid Helium-32&2£. 
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For a two-dimensional rotationally invariant system 
with quadratic dispersion, 



m 

X{<1) = 2~, 



1 - 



ReyV - (2k f y 



(16) 



Note that the susceptibility is a constant for q < 2k f. 
However, for a system with a circular Fermi surface, the 
furthest apart any two points on the Fermi surface can 
be is 2k f. Therefore, the pairing strengths correspond 
to the eigenvalues of a constant matrix. One sees imme- 
diately from this that the s-wave solution has a positive 
eigenvalue A = m/27r+ci, and all higher angular momen- 
tum gap functions have zero eigenvalues. Therefore, for 
such a rotationally invariant system in two dimensions, 
superconductivity does not occur in any channel to order 
U 2 . When higher order terms in perturbation theory are 
taken into account, it has been found that the leading 
instability occurs in the p-wave channel^. 




FIG. 1: Density of states as a function of electron concentra- 
tion on the square lattice for t' = (solid line) and t' = —0.3 
(dashed line). 



B. Lattice systems: d=2 

Various authors^— have studied the Hubbard model 
on the 2D square lattice in the weak coupling limit. Near 
half- filling, nesting effects lead to d x 2_ y 2 pairing. In weak 
coupling at low doping, higher order non-quadratic terms 
in the quasiparticle dispersion^ lead to d xy pairing for 
t' = and for t' < —t/A, the p-wave state is favored^. 
In higher order perturbation theory^, U 3 vertex correc- 
tions have been shown to enhance the p-wave over the 
d X y state. Here, we are interested in the behavior of the 
model as U jt approaches zero while keeping the density 
of electrons fixed, and we have studied the eigenvalue 
problems given in Eqs. [7] and [5] must be studied numer- 
ically. A finite number of points on the Fermi surface 
form the basis for the matrices g s ' 1 . These in turn are 
diagonalized and the pairing eigenvalues are determined 
for each electron concentration. For each of the systems 
considered below, we have found that a discretization of 
the Fermi surface with 500 — 800 points is more than suf- 
ficient to produce accurate results. Both spin and the 
lattice point group symmetries are used to classify the 
resulting eigenstates. For a tetragonal crystal (i.e. the 
point group D^) without spin-orbit coupling, there are 
4 non-degenerate spin singlet states, and a two-fold de- 
generate triplet state: 
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where the left hand side labels the irreducible represen- 
tation and the right hand side lists the basis functions 
(with the association x — > sirik x ,x 2 — > cosk x , etc). 
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FIG. 2: Pairing strengths for the 2D Hubbard model at t' = 
as a function of electron concentration. 

Figure Q] shows p as a function of electron concentra- 
tion n for a tight-binding model on a square lattice. Fig- 
ure [2] shows the pairing strengths for the 2D Hubbard 
model on a square lattice with t' = as a function of n. 
(Particle-hole symmetry assures that the phase diagram 
is invariant under n — > 2 — n.) We see clearly that near 
half- filling {n = 1), the dominant form of superconductiv- 
ity has d x 2_ y 2 symmetry. At about n = 0.6, we see that 
the favored configuration changes to d xy pairing. Thus, 
there are two distinct superconducting ground states that 
occur on a square lattice with a near-neighbor hopping 
as a function of concentration for the particle-hole in- 
variant system at t' = 0: there is a d x 2_ y 2 ground state 
for 1 > n > 0.6, and a d xy ground state for n < 0.6. 
The results shown in Fig. [2] are in qualitative agreement 
with the spin fluctuation exchange studies of the Hub- 
bard model of Scalapino et a$&>. Figure [3] shows the gap 
function for the d-wavc state at n — 1.1. Notice that 
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FIG. 3: The gap function A for the d x 2_ y 2 state which occurs 
at t = and n — 1.1. The pair field is plotted as a function 
of 9f, the angular degree of freedom on the Fermi surface 
relative to the x axis. The pair field obtained in the weak 
coupling analysis deviates significantly from a simple cos(2#p) 
form, depicted here by the dotted line. 




FIG. 4: Pairing strengths for the 2D Hubbard model at t' = 
— 0.3i as a function of electron density. The critical density 
n c at which the van Hove singularity occurs at the Fermi level 
is shown. 



it has a shape that differs substantially from the simple 
(cosk x — cos/cj,) form. 

For t' 7^ 0, on the square lattice, the particle- hole 
symmetry is destroyed and the van Hove singularities 
occur away from half-filling. Figure H] shows the pair- 
ing strengths on the square lattice at t 1 = —0.3 as a 
function of doping. Again, the dominant configuration 
which occurs near the half- filled system is d x i_ y i pair- 
ing. The van Hove singularity occurs in this system at 
n = n v h « 0.72. In a narrow window of densities near 
n V H, the d-wave order is suppressed and the p-wave order 
dominates. Upon further decreasing the electron concen- 
tration away from the van Hove singularity, the d x 2_ y 2 
ground state again gives way to a p-wave superconduct- 
ing state at n k 0.5. Qualitatively similar results have 
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FIG. 5: The effective interaction for t' — —0.3 (dashed and 
dashed-dotted lines) and for t' — —0.05 (solid line) as a 
function of n — At' . For t' = —0.3, the dashed line repre- 
sents the pairing strength of the d x 2_ y 2 state whereas the 
dashed-dotted line corresponds to the pairing strength of the 
p-wave state. While both are depressed as the chemical po- 
tential crosses the van Hove singularity, the p-wave state ob- 
tains a larger pairing strength parametrically close to the van 
Hove point. For t' — —0.05, the nesting of the Fermi sur- 
face remains nearly perfect and acts to enhance the pairing 
strength of the d x 2_ y 2 state. The p-wave pairing strength for 
t' = —0.05 is not shown here, since it is much lower than the 
d-wave strength. 



been found for the square lattice t, t' Hubbard model us- 
ing functional renormalization group analysis for fillings 
close to the van Hove singularities^. 

It is worth examining the singular behavior near n v n 
(i.e. near the point /j = At') in more detail. For the 
square lattice with t' ^ 0, the susceptibility at q = (0, 0) 
diverges logarithmically as ji — > At 1 : 
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For a large momentum transfer, Q = (ir,ir), the suscep- 
tibility varies as 
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for |/i — At'\ << t'. Thus, a finite t' acts to sup- 
press the nesting of the Fermi surface ( for the perfectly 
nested Fermi surface at t' = 0, x{Q) diverges much more 
strongly as In 2 Taking the dominant scattering 

processes at momenta q = (0,0) and Q = (tt,tt) into ac- 
count, one finds that the effective d-wave pairing strength 
near the van Hove singularity is roughly 
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where Vq is a subdominant contribution which arises from 
the intermediate momentum transfers on the Fermi sur- 
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FIG. 6: (upper left): Energy contours of the triangular lat- 
tice nearest-neighbor tight-binding model. The dashed black 
line marks the zone boundary. Blue contours correspond to 
the band bottom and red contours occur near the top of the 
band, (upper right) Density of states on the triangular lat- 
tice. When the Fermi level is at the van Hove point, the 
volume of the Fermi surface is 3/4 of the zone, (lower left) 
Energy contours of the honeycomb lattice nearest-neighbor 
tight-binding model. Due to the particle-hole symmetry, only 
the e > contours are shown, (lower right) Density of states 
of the honeycomb model. The van Hove singularities occur at 
a filling of 3/8 and 5/8. 




FIG. 7: The two distinct f-wave irreducible representations 
on the triangular and honeycomb lattices, (left) Bi u : ip ~ 
k y (fcy - Ski), (right) B 2u : ip ~ k x (k% - 3fcJJ). The Bi u gap 
function is the dominant gap when there are two disconnected 
Fermi pockets in the system, each centered around the zone 
corners. The Bin gap function has the opposite sign on each 
pocket. 



face. For t' << t, the d-wave pairing strength is en- 
hanced as \fi — 4t'| — > 0. However, for t/2t' > e the 
d-wave pairing strength decreases as the van Hove singu- 
larity is approached. This behavior is illustrated in Fig. 
[5l It reflects the fact that the d-wave pairing is driven 
by an interaction which is stronger at large momentum 
transfer. An interaction which is greater at small mo- 
mentum transfer suppresses the d-wave pairing. 

Next we consider the nearest-neighbor tight-binding 
model on the 2D triangular and honeycomb lattices. Fig- 
ure [5] shows the basic electronic structure on these lat- 
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FIG. 8: The dominant pairing strengths on the triangu- 
lar lattice as a function of electron concentration. The two- 
component d-wave order parameter occurs in the hole-doped 
(i.e. x < 0) system. But for electron-doping (x > 0) beyond 
the van-Hove filling, the f-wave(2) gap develops the largest 
critical temperature. 

tices. Both lattice systems have the hexagonal point 
group (Deh ) symmetry and therefore the same irreducible 
representations characterize their gap functions. In the 
singlet (even parity) channel, the following are the irre- 
ducible representations: 

Singlet channel 

A lg : ?P~k 2 x + k 2 y 

E 2g : i>~{(kl-kl),{2k x k y )} 
A-2 g ■ ip ~ k x k y (k 2 — 3ky)(ky — 3k 2 ) 

(21) 

Note that the d-wave function is a two-dimensional rep- 
resentation (this is in fact the only two dimensional rep- 
resentation in the singlet channel for this system, so long 
as we restrict our superconductivity to be only in the 
basal plane). For the triplet channel, we have 

Triplet channel 

Eiu ■ ip ~ {k x , k y } 
B lu : ip ~ k y (ky- ik 2 x ) 
B 2u : i> ~ k x (k 2 x - 3k 2 y ) 

(22) 

The p-wave gaps form a two-dimensional irreducible rep- 
resentation. There are two distinct one-dimensional f- 
wave gap functions that belong to the B itl , Ji 2 u repre- 
sentations and are shown in Fig. [7] 

Figs H] and M show the pairing strengths for the tri- 
angular and honeycomb lattice respectively. For both 
systems it is seen that the dominant pairing instabilities 
are either the two-component d-wave representation, or 
the non-degenerate f-wave representation. 
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FIG. 9: Pairing strengths for the honeycomb lattice. The 
trends on the honeycomb lattice are some what reverse in 
comparison to the triangular lattice. However this can be 
understood from looking at the Fermi energy contours. For a 
given sign of t, the energy contours of honeycomb lattice for 
large electron doping (i.e. top of the band) is similar to those 
of the triangular lattice at large hole doping (i.e. bottom of 
the band) and vice-versa. 



For the triangular lattice at modest electron concentra- 
tion, or the honeycomb lattice far away from half-filling, 
the Fermi surface is hexagonal in shape and is simply con- 
nected. In this regime, the dominant pairing configura- 
tion is d-wave pairing that consists of a d x i_ y i component 
degenerate with a d xy component^ In our analysis, the 
magnitude of each of the two d-wave gap functions and 
their relative phase cannot be determined, since they are 
determined by non-linear effects captured, for instance in 
the Landau Ginzburg theory. However, it is reasonable 
to expect that in order to gain condensation energy, the 
system will spontaneously break time-reversal symmetry 
and form a d + id superconductor. 

For the triangular lattice near the top of the band, 
and for the honeycomb lattice close to half-filling, the 
Fermi surfaces form disjoint pockets. In this concentra- 
tion regime, we have found for both systems that the 
triplet f-wave gap function is the ground state. The f- 
wave state which is favored has its lines of nodes along 
the lines connecting the zone center to the midpoints of 
the zone edges (see Fig. [7]). However, since these lines 
of nodes never cross the Fermi surfaces which are cen- 
tered on the corners of the Brillouin zone, the f-wave 
pairing produces a fully gapped state on the Fermi sur- 
face. The gap changes sign between the two distinct 
Fermi surfaces. It has been argued in the past based on 
the spin-fluctuation exchange mechanism that the strong 
magnetic excitations associated with such disjoint, rea- 
sonably well-nested Fermi pockets results in an effective 
pairing interaction which is repulsive between the two 
pockets. Consequently, a solution which encodes a sign 
change of the gap among the two Fermi surfaces will nat- 
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FIG. 10: Density of states as a function of electron concentra- 
tion for the diamond lattice (solid black curve), SC (dashed- 
dotted curve), BCC (dashed curve) and FCC (blue dashed 
curve) lattices. The FCC is a non-bipartite lattice and there- 
fore, the particle-hole symmetry is absent. 

urally be favored^. The transition between the f-wave 
and d-wave pairing states occurs as the electron density 
is varied across the van Hove filling. We note, finally, 
that the f-wave solution has a substantially higher tran- 
sition temperature than the other superconducting gap 
functions found in these lattice systems. 

C. Lattice systems: d=3 

Next we consider 3 dimensional lattice systems. In 
particular, we shall consider the superconducting insta- 
bilities of the Hubbard model on the simple cubic (SC), 
body-centered cubic (BCC), face-centered cubic (FCC) 
and diamond lattices. We shall restrict our analysis to 
nearest neighbor tight binding dispersion in each of these 
cases. All of these lattices have the octahedral (Oh) point 
group symmetry which permits the following classifica- 
tion for the gap functions: 
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(23) 

for the singlet gap functions and 

A Xu : ip ~ xyz 
Tin ■ tp ~ {x, y, z} 

T 2u :V ~ {x(y 2 -z 2 ),y(z 2 -x 2 ),z(x 2 -y 2 )} 

(24) 

for the triplet states. Figure [TU] shows the density of 
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FIG. 11: Pairing strengths of the dominant pairing configura- 
tions on the simple cubic lattice. For electron concentrations 
that are far from half-filling, the 3-fold degenerate p-wave 
(Tiu) state occurs whereas closer to half-filling the doubly- 
degenerate d-wave state (E 9 ) is found. 



0.08 







— {xy,yz,zx} 




-{x 2 -y 2 ,3z 2 -r 2 } 




— {x,y,z} 

\ 

** ♦** * 




• . . • " " 


» f $ 
***♦ 







0.2 0.4 0.6 0.8 

n 



FIG. 12: Pairing strengths for the dominant pairing configu- 
rations on the BCC lattice. For electron concentrations that 
are far from half-filling, the 3-fold degenerate p-wave (Tiu) 
state has the greatest strength. For intermediate concentra- 
tions, the E 9 and T2 9 d-wave solutions are found to have the 
strongest strengths. 



states as a function of electron concentration on each of 
these lattice systems. With the exception of the BCC and 
FCC lattices, the density of states remains finite across 
van Hove singularities in 3 dimensional systems. 

Figures [TT] and [12] display the phase diagram on the 
SC and BCC lattice respectively as a function of elec- 
tron concentration. At low concentrations, the p-wave 
solution has the highest T c in both lattices whereas near 
half-filling, a d-wave solution has the higher transition 
temperature. On the SC lattice the E ff d-wave configu- 
ration is found near half filling. This state is the 3 dimen- 
sional analog of the d x 2_ y 2 gap found on the square lat- 



FIG. 13: Pairing strengths for electrons on a FCC lattice. 
As in the other 3D lattice systems, p-wave superconductiv- 
ity occurs for low electron concentrations. For intermediate 
concentrations, the d-wave solutions have the largest pairing 
strengths. Near the top of the band, the density of states 
diverges and a 3-component f-wave superconductor has the 
strongest pairing strength. 



tice near half-filling. Due to the underlying symmetry of 
the cubic lattice, the d x 2_ y 2 gap must be degenerate with 
the d z 2_ r 2 gap function. This degeneracy in turn implies 
that below T c , the system on the cubic lattice will spon- 
taneously break time- reversal symmetry and form a d+ id 
gap function. On the BCC lattice (Fig. [T2"]) . the dom- 
inant pairing configuration near half-filling is the triply 
degenerate d-wave gap function (T 2g ). At intermediate 
concentrations, the E g gap functions have the largest T c 
and again at low concentrations, the p-wave solution is 
favored. On the FCC lattice (Fig. the T 2u f-wavc 
gap is favored near the top of the band. This f-wave state 
gives way to the E g d-wave state below n = 1.5, which 
in turn is replaced by the T 2g d-wave configuration be- 
low n = 1.3. This triply degenerate d-wave gap function 
persists for a wide range of concentrations. Finally, at a 
concentration below n = 0.1 we find the p-wave state has 
the highest pairing strength. Finally, on the diamond lat- 
tice (Fig. [14]), the T 2u gap occurs near half filling when 
the semi-metal is lightly doped. As the concentration is 
increased, the ground state consists first of the E ff and 
subsequently the T 2g superconducting states. Near the 
top and bottom of the band, the p-wave solution has the 
largest pairing strength. 



IV. RG STRATEGY 

At any given temperature, the thermodynamic prop- 
erties of the model can be computed perturbatively in 
powers of U /t so long as U/t is small compared to a 
characteristic T dependent magnitude. (Indeed, we gen- 
erally expect that perturbation theory is convergent at 
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FIG. 14: Pairing strengths for electrons on a diamond lat- 
tice. For low electron concentrations, the p-wave solution 
has the largest pairing strength whereas near half-filling, the 
triply degenerate f-wave gap has the greatest pairing strength. 
For intermediate concentrations, both the doubly and triply 
degenerate d-wave pairing states have the dominant pairing 
strengths. 



finite temperature, with a finite radius of convergence. 
Dynamical properties of the system often depend non- 
analytically on the strength of the interactions, even at 
elevated temperatures, but as long as we stick to ther- 
modynamic quantities, this issue should not arise.) Al- 
ternatively, even at T = 0, if we introduce an artificial 
low energy cutoff, fio, m the spectrum, low order per- 
turbation theory is reliable so long as fl > Qpt, where 
ilpT can be obtained by looking at the most divergent 
terms in each order of perturbation theory, the familiar 
particle-particle ladders: 



p\U\\og[W/n PT } = 1; 

n PT = W exp[-l /p{E F )\U\]. 



(25) 



ClpT is a physical energy scale in the problem - it is the 
highest energy at which the bare interactions begin to be 
significantly renormalized by many-body effects. 

With this in mind, we formulate the problem in terms 
of a Grassman path integral, which we express in terms 
of the normal modes of the quadratic piece of the action 
defined by Ho. As a first step, we integrate out all the 
modes with energies greater than a cutoff, chosen so 
that 



U 2 /t > fi > SIpt- 



(26) 



Because Qq 3> &>pt, the interactions in the resulting ef- 
fective action can be computed using straightforward per- 
turbation theory. Moreover, we are guaranteed that in 
dimensionless units, all the effective interactions are still 
weak. Because W U 2 /t, the resulting effective action 
involves only modes within a paramctrically narrow win- 
dow, of width Qq, about the Fermi surface. In particular, 



this effective action is of precisely the form assumed as 
the starting point for the perturbative RG analysis of 
the Fermi liquidi 31 ' 32 Specifically, the dispersion can be 
linearized about the Fermi surface, the effects of small 
irrelevant terms can be neglected, and the beta function 
for the marginally relevant interactions can be computed 
to one loop order. 

The second step is to compute the RG flows starting 
from the initial data obtained in the first perturbative 
step. These flows describe how the effective couplings 
change as we continue the process of integrating out high 
energy modes by reducing the cutoff below SIq. These 
equations cease to be accurately governed by the per- 
turbative beta function when one or more dimensionless 
interaction grows to be of order 1. However, the value of 
the cutoff, ft*, at this point defines (up to a multiplicative 
constant of order 1), a characteristic energy scale in the 
problem. Assuming that a one parameter scaling theory 
describes the low energy physics, then all emergent ener- 
gies in the problem, including T Cl the root-mean-squared 
gap magnitude, Ao, etc., are all simply proportional to 
£1*. (Without knowing more explicitly the crossover be- 
havior from the Fermi liquid to the superconducting fixed 
point, it is not possible to obtain a precise value for the 
constants of order 1, and hence ao in Eq. [5] cannot be 
computed by the present methods.) 

Note that Oo in this treatment is not a physical en- 
ergy scale, but rather a calculational convenience. It is 
important that the results should be independent of the 
value of £lo- We will see that by chosing f^o <C U 2 /t, 
we make simple the manipulations that insure that our 
results are independent of f^o, at least to the desired or- 
der in powers of U/t. This highlights an important dif- 
ference between the present analysis and a conceptually 
similar treatment of the small U problem considered by 
Zanchi and Schula^, in which a two step RG analysis 
of the Hubbard model was undertaken, but in which the 
intermediate scale (which they call £ x ) is a physical U in- 
dependent measure of the proximity to a van-Hove point. 
(We return to this comparison in Sec. IVII1 below.) 



V. FIRST STAGE RENORMALIZATION: 
PERTURBATIVE RESULTS 

The first step is to integrate out the states with ener- 
gies down to fio, and to compute the resulting effective 
interactions perturbatively in powers of U. The impor- 
tant terms, which will serve as inputs for the second stage 
renormalization, are the electron self energy, £(£;), and 
the two particle vertex, T a>a > (fc, q), in the particle-particle 
channel i.e. the scattering amplitude of a pair of particles 
with spin polarization a and a' and momenta k and — k 
into a pair of particles with the same spin polarization 
and momenta q and — q. The two particle vertex can be 
decomposed into the singlet and triplet channels, r s and 
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r 4 respectively, where 



-t t- 
+ 

-I I- 



-t t- 
+ 

-I I- 



T s (k,q) = - [T^(k,q)+T U (k,q) 



W T t (k,q) = r t , t (fc,g) (27) 

( 2a ) ( 3a ) ( 3b ) The electron self-energy and two-particle vertices are 

t ' t t ' g » t t > ; » t tWt 



(28) 



(3c) 



(3d) 



(3e) 



(3f) 



t 



T-^t t— r— t 

6+ + A+ + 
f-» — ^-»-f f-*-^ — ^»-t \t 



(2b) 



(3g) 



(3h) 



FIG. 15: Two particle irreducible diagrams to third order in 
U which contribute to the effective interaction in the particle- 
particle channel. The incoming set of electrons at the left of 
each diagram have momentum k, —k and are scattered by 
the interaction to states with momenta q, —q. The dashed 
line denotes the bare Hubbard interaction, and the solid lines 
correspond to the bare electron propagators. The diagrams 
are constrained by the requirement that the spin must be 
flipped across the dashed line. 
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FIG. 16: Diagrams to 0{U 4 ) which are used in determining 
the effective interaction in the particle-particle channel. The 
upper set of diagrams (L a -Lfc) show the contributions in the 
spin singlet channel, whereas the lower set (L m , L n ) show the 
processes that contribute in the triplet channel. 



s (2) (fc;^c 



u 2 ^ 2 (k;n Q ) + o{u 3 ) 

^L-G(k + q^)G(q^) (29) 
fin ( 2lT ) ' 



r s (k,q-,n ) = 

r t (k,q;tt ) = 

rf\k,q;Q ) = 

P(O ) = 



U + Y^U n T^(k,q;Q ) 

n>2 



(30) 



(31) 



Tl>2 



x(k + q;n ) + P(Q Q ) (32) 
-x(k-q;Q ) (33) 

i^G(p>)G(-p,-u/)(34) 
o (27r) d+i 



Here, G(k; w) is the single-particle Green function, J" n 
signifies the integral over all q subject to the constraints 
le^l > J7o, and > ^o, and d is the number of spa- 

tial dimensions. The first term in Eq. [35] is obtained 
from diagram (2a) of Fig. HT)1 whereas the second term is 
obtained from diagram (Lf,) in Fig. 1161 The quantity in 
Eq. is obtained from diagram (2b) in Fig. QjjJ Below, 
we shall discuss the higher order terms, T^™' with n > 2. 

The particle-hole bubble, x, is regular in the limit 
Slo — > 0, while the particle- particle bubble has the 
well known logarithmic divergence associated with the 
Cooper instability, but is otherwise regular: 



X(fc;0o)= X (k) + O(Cl ) 

p(n )= P iog[A/n ] + o(n ) 



(35) 



where A is determined by the band-structure over the en- 
tire band. Since we will never need to keep terms higher 
order than U 4 , and since, by assumption, il <C U 2 , 
the higher order terms in powers of f2o can be neglected 
henceforth. 

The higher order vertex functions can likewise be eval- 
uated by keeping a non-zero value of f2o where ever there 
is a logarithmic divergence in the fio — ► limit, but set- 
ting f2o = elsewhere. The expression in the singlet 
channel is 



if )(fc,(f;Oo) = p 2 log 2 [A/n ] + [ 7 < 3 >(fc) + 1 ^(q)]plog[A/n ] + f?\k,q) + O(fio) 



(36) 
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where if we define p to designate a vector on the (un- 
perturbed) Fermi surface and Sf = J dpp equal to the 
"area" of the Fermi surface 

^® m f *,(■%>) *<! +f> (37) 

with vp (p) the magnitude of the Fermi velocity at posi- 
tion p on the Fermi surface and the norm of the Fermi 
velocity defined according to 

^../*(» y (38) 

vf J of \v F {p) J 

The first term of Eq. [35J is obtained from diagram (L^) 
in Fig. [TBI whereas the terms proportional to log [A/f2o] 
are derived from diagrams (L e ,L/) in Fig. [15] with the 



thick solid line treated to second order in U. Finally, f ( 3 ' 
contains all the non-singular contributions in the limit 
f2o — >• 0, which are derived from the third-order two- 
particle- irreducible (2pl) diagrams (3a-3f in Fig. [15]) . 
They can be expressed as double momentum integrals 
over suitable products of quartets of G's, as shown in the 
appendix, but in the interest of clarity, we do not display 
them here. In the triplet channel, the third-order correc- 
tion to the vertex is non-singular and is obtained from 
diagram (3g) of Fig. [T5j 

T?\k, q; fi ) = f [ 3 \k, q) + O(fio) (39) 

An explicit expression for this quantity is also derived in 
the appendix. 

Similarly, we can obtain an expression for r*- 4 ^ 



ri 4 >(fc,q;tt ) = {p 3 log 3 [A/n } + [ 7 < 4) +1 [ 3 \k) + 1 [ 3 \q)]p 2 log 2 [A/Q ] 

+ [l i 2 l} (k) + 7 ^ 4) (q1]plog[A/O ] + 7 i 4) (fc, q) P \0g[A/n Q }} + f^(k, q) + O(Q ) 

I 



(40) 



where 



(4) 



dpdp' 



Si v f (p)vf(p') 



^(k,q) 



F 

dp f vp 



fW(k,p) 



S F \Vf(P) 

1Tf x{Up) (^§) 



(41) 
(42) 



X (q+P) (43) 



The first term in Eq. 05] is represented by the fourth 
order ladder diagram in (L g ) of Fig. 1161 The terms 

involving 7^ are obtained from diagrams (L/j-Ly), and 

those involving 7^ 4 ' are derived from diagrams (L e ,L/), 
treating the thick solid line to 3rd order in U. Finally, 
f i 4 -* is a non-singular contribution from the fourth-order 
2PI diagrams, which are not shown here Since we shall 
not make use of these terms, we will not provide explicit 
expressions for them. In the triplet channel, the fourth- 
order vertex consists of a single term 

r| 4) (fc,q;r! ) = ploglA/iloh^ik^ +O(n ) (44) 
where 



VF 
Vp{p) 



xip- q) (45) 



which is obtained from diagram (L„) of Fig. [T5J treating 
the thick solid line to second order in U. 



VI. SECOND STAGE RESULTS: RG ANALYSIS 

The second stage of renormalization is carried out fol- 
lowing the standard Fermi liquid renormalization group 
(RG) procedure of Shankar— and Polchinski^. Notice 
that the effective action generated after the first stage of 
renormalization is precisely of the form assumed as the 
starting point of this renormalization procedure: the ef- 
fective interactions are all small (so a perturbative RG 
approach is justified) and the remaining states lie in a 
narrow strip of width fio about the Fermi surface so 
that the spectrum can be linearized without loss of ac- 
curacy. Various interactions (Fermi liquid parameters) 
are marginal at the non-interacting fixed point. These 
do not significantly affect our principle results. The only 
couplings that renormalize are those in the Cooper chan- 
nel. These are governed by the one-loop RG equations 
which we write in matrix form as 



dg 

di 



= -.9*5 



where I = log[f2o/fi], 



dp , 



where g designates a dimensionless matrix 



9i, 



i vf 

VF{k) 



T(k,q) 



VF 
VF(q)' 



(46) 



(47) 



(48) 



Here, we have left implicit the dependence of both g and 
r on the spin indices. 
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r 



(b) 



FIG. 17: An arbitrary irreducible interaction vertex T(k, q) in 
the particle-particle channel is shown in (a). The spin indices 
are suppressed for clarity. The flow of the most negative eigen- 
value of F breaks down at a scale which is identified as the su- 
perconducting transition temperature. For each such vertex, 
there belongs a corresponding diagram of the form shown in 
(b), which produces the logarithmic divergence ~ Iog(j4/Oo) 
needed to remove the dependence on the arbitrarily chosen 
initial cutoff f2o. 



In integrating these equations, we start with an initial 
value of the interaction matrix, g°, which is the output 
of the first stage of renormalization. Because g is a real 
symmetric matrix, it is also Hermetian, and so can be 
diagonalized: 



4 ? =E A °^(W»(<7) (49) 



where the eigenvalues, A„, are real and the cigcnfunctions 
form an orthonomral basis, 



dk 



CO)V'm(fc) = Sn 



(50) 



As a result, each eigenvalue renormalizes independently 

A„ 



(51) 



To determine the physically important scale, f2*, at 
which the RG treatment breaks down, we must first iden- 
tify the smallest eigenvalue of g°, Aq 1 , for which Ag < A° 
for all n > 0. Assuming that Aq is negative, then 



n* = n exp[-l/|Ag| 



(52) 



Asymptotically, both T c and the zero temperature gap 
scale, Ao, are then equal to a (unknown) number of or- 
der 1 times 0*. Under generic circumstances, we expect 
Ao to be non-degnerate, i.e. Xq < X\. The exception to 
this is the vicinity of a zero temperature phase transition 
between two superconducting states with different sym- 
metries, where two eigenvalues cross. The properties of 
the infinitesimal neighborhood of such critical points will 
not be investigated further in this paper. 

There is an apparent problem with Eq. 1521 which is 
that it has an explicit dependence on fio- Since fig was 
introduced as an unphysical calculational device, this de- 
pendence must be spurious. Fortunately, there is also an 
implicit dependence of Aq on fl^, which just cancels this 
explicit dependence. 

To see this most simply, first consider the case of the 
negative U Hubbard model. In this case, the lowest 



cigncvalue is clearly in the spin singlet channel, and it 
can be computed perturbatively in powers of U as 



l v F (k) 



[l + 0(£/)] 



x u — 



V vf 

P U + pU 2 [P(n )+p'} + O(U 3 ) (53) 



where 



P' = 



dkdq . - „. 
-52-X(k + q). 

iJ TP 



(54) 



The logarithmic dependence of P on f2rj insures that 
when the expression in Eq. l5Hlis inserted into Eq. 152] the 
result is independent of SIq, at least to the stated order 
in powers of U: 



Aexp[-(l/p\U\)e-"' /p \l + 0(17)1 . (55) 



The same analysis can be carried through for the more 
complicated case of the repuslive U Hubbard model. In 
this case, however, since the leading contribution to Aq 
is order U 2 , the term responsible for canceling the JIq 
dependence of the prefactor is a logarithmically divergent 
fourth-order order contribution to the 'dressed' vertex. 
The dressed vertex is represented by the thick solid line 
in diagram (L c ) or (L m ) of Fig. [TBI We shall consider the 
singlet and triplet channels separately in what follows. 

In order to see that fT in the singlet channel does not 
depend on Slo, it is helpful to discrctizc the points on 
the Fermi surface, so that the matrix, T s in Eq. 1321 and 
hence the matrix, g s in Eq. 551 as well, is a TV x TV matrix, 
where TV is the number of k-points. (The continuum limit 
can easily be taken at the end of the calculation.) This 
TV x N matrix can then be partitioned into a 1 x 1 block 
g Si o which affects only the trivial s-wave states, and an 
(TV — 1) x (N — 1) block <? S: i containing the non-trivial 
superconducting states: 



g s {ki, kj) 



g s .o T s 
V 9s.i 



(56) 



Here, T s is a 1 x (TV — 1) matrix which connects the 
trivial s-wave subspace with the orthogonal space of sign- 
changing pair-fields. 

The lowest order contributions to g s ,o oc r Si o comes 
from diagrams (L a ,Lf,) in Fig. [T6jso T s ,q ~ 0(U). On the 
othcrhand, g Sy \ oc T Sf i ~ 0(U 2 ). The diagrams which 
contribute to Ts(ki,kj) are shown in (L e ) and (L^) in 
Fig. [TH1 The interaction vertex in these diagrams has 
the form of a product of two terms, one of which consists 
of a bare (and therefore momentum-independent) vertex 
operating on the incoming momenta (hi), and the other 
of which is a dressed vertex connected to the outgoing 
momenta (kj). Therefore, these terms, when viewed as a 
matrix, operate on a non-trivial singlet pairing configura- 
tion and produce the trivial s-wave solution. From Figs. 
[TSIandfTfjl it is easy to see that T s ~ 0(U 3 ). In a similar 
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fashion, diagrams (L/) and (Lj) contribute to 7~J which 
is just the hermitian conjugate of T s . In diagonalizing 
g s , the off-diagonal matrix elements, T s , can be treated 
perturbatively. Indeed, it is clear that the leading effect 
of these terms on the eigenvalue problem is 0(U 5 ), and 
so to the order we are working, they can be set to zero. 

Within the non-trivial (N — l)-dimcnsional subspace, 
the lowest-order term in perturbation theory which con- 
tribute to r Sj i is represented by the dressed vertex in 
diagram (L c ) of Fig. [T6l To lowest order in U/t, this is 
the diagram in (2a), which is U 2 x(ki + kj ;; ^o)j an d when 
all of diagrams (2a)-(3g) are taken into account, the non- 
singular terms of 0(U 3 ) which produce fi (hi, kj) in Eq. 
(3B]will also contribute to the vertex. We define the quan- 
tity T s i to be the piece of r Sj i which is non-singular in 
the limit i7o — > 

f sA (k u kj) = U^ik+k^ + U^Kkk^+OiU 4 ) (57) 

and correspondingly 

9s,i(k, q) = p\Jv F /v{k) f s ,i(fc, q) y/v F /v(q). (58) 

The lowest order term in perturbation theory for T 
which has singular (logarithmic) Hq dependence is the 
same diagram (L b ) in Fig. [T5] that gives the logarithm 
in Eq. 1531 however, this term is purely a contribution 
to r s .o- The lowest order singular term in r s> i which 
has the singular ft dependence derives from fourth-order 
diagram (L&) in Fig. [TH Thus, when all the diagrams 
are properly taken into account, we see that 

9s,i = g.,i + g a ,i * g.,i log [A/fio] + 0(U 5 ) (59) 

Thus, if we express the results in terms of A, the eigen- 
values of the non-singular part of the interaction, g s i, we 
find that 

T c = Acxp{-l/|A mm |} [l + 0(U 2 )] (60) 

The logarithmic dependence on ilo of the fourth order 
contribution to Aq has just the requisite form to cancel 
the explicit dependence on Qq- 

The eigenvalue obtained this way is valid to 0(U 3 ). 
However, to obtain the quantity a^, one only needs to 
consider g s ^ to 0(U 2 ), which is obtained from the first 
term in Eq. 1571 Once the eigenvalue problem to this or- 
der has been solved, the 0(U 3 ) correction, which we refer 
to as a.\, is obtained by treating the contribution from 
the second term in Eq. [57] to first order in perturbation 
theory. 

It is similarly straightforward to show that T c for a 
spin-triplet ground state is independent of the initial cut- 
off Qq . Again, the effective interaction in the triplet chan- 
nel must have non-trivial momentum dependence, and 
only the dressed vertex shown in Fig. [T5]can contribute. 
The resulting expression in the triplet channel is directly 
analogous to Eq. [57J 

f t.iih h) = -u 2 x(k - h) + u 3 r?\ki, h) + o(u 4 ) 

(61) 



The dependence of T c on the initial scale fio is elimi- 
nated by diagram (L„) in Fig. [TBI which possesses the 
required logarithmic divergence ~ log[A/fio]- After in- 
cluding both diagrams (L m ) and (L n ) in Fig. 1161 we find 
that the interaction vertex in the triplet channel is 

gt,i = hi -gt,i* h,t log [A/n ] + o(u 5 ) (62) 

Thus, as was found in the singlet channel, the final ex- 
pression for T c is independent of Qq. 

The general feature in perturbation theory which acts 
to eliminate the initial cutoff dependence of T c is shown 
in Fig. 1171 In Fig. 1171 (a). an arbitrary irreducible inter- 
action vertex T in the particle-particle channel is shown. 
For the negative U Hubbard model, T would correspond 
simply to the bare vertex, and for the repulsive cases, 
r is the appropriate irreducible vertex in either the sin- 
glet or triplet channel as discussed above, and could be 
computed to an arbitrary order in perturbation theory. 
For each such V, there corresponds a diagram of the form 
shown in !17f b) which acts to remove the dependence on 
f^o- In this diagram, the internal legs which separate 
each vertex T produces the required logarithmic diver- 
gence factor P(£lo). Since this is true for an arbitrary 
interaction vertex, the fact that the final expression for 
T c is independent of the arbitrarily chosen initial scale 
flo is true to all orders in perturbation theory for the in- 
teraction vertex (we note in passing that self energy cor- 
rections to the internal propagator lines produces 0(U 2 ) 
corrections . in the Hubbard model and do not affect the 
general structure discussed here). 



VII. RELATION TO PREVIOUS WORK 

The present work is, in part, a recasting of old work 
in a new framework, in terms of a more well controlled 
asymptotic analsyis. Specifically, Kohn and Luttinger— 
observed that even for a repulsive bare interaction, 
the momentum-dependent structure in the irreducible 
particle-particle vertex can give rise to a Cooper insta- 
bility in a suitable channel. For a short range repulsive 
bare interaction, the resulting pairing interaction is me- 
diated by an S = 1 particle-hole channel which is the 
leading term of the Berk-Schrieffcr spin fluctuation ex- 
change interaction^. 

An important piece of the physics of correlated ma- 
terials that is missing in the weak coupling limit is as- 
sociated with "competing orders," and the accompany- 
ing interplay between interactions in different "channels." 
We stress that this is not a failure of the method of 
solution, but is something that is an intrinsic feature 
of a weakly perturbed Fermi liquid. In order to ad- 
dress this physics, various calculations have been un- 
dertaken using the Functional Renormalization Group 
(FRG) metho d 17 ' 25 " — which is in many ways similar in 
structure to the weak-coupling analysis undertaken here. 
In this approach, a single set of pcrturbative non-linear 
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flow equations are derived for the coupling constants de- 
fined both near to and far from the Fermi surface, and 
then the resulting flows are analyzed numerically starting 
from initial conditions corresponding, for instance, to the 
Hubbard model with intermediate couplings, U ~ 3 — At. 
As stressed in Ref^, for these values, the effective in- 
teractions near the Fermi surface are typically large. On 
an intuitive level, the advantage of the FRG approach is 
that it does capture some physics of multiple intertwined 
scattering processes in a physically compelling manner. 
On the other hand, the perturbative methods used are 
formally valid only in the asymptotic limit U/t — > 0. 
While, as far as we know, there has been no published 
work analyzing the FRG flows for the Hubbard model in 
this limit, were the FRG calculations carried out in this 
limit they would approximate the leading order behav- 
ior derived in the present paper using somewhat different 
methods. 

Another feature of the weak coupling limit is that the 
dynamics is set by the bare bandwidth. Thus, the bare 
susceptibilities that enter the expression for the pairing 
interaction are those of the unperturbed Fermi gas, in 
which the only energy scale is the bandwidth. Again, a 
plausible extension of the present results involves replac- 
ing one or more of these susceptibilities with a dressed 
susceptibility, or even with an experimentally measured 
susceptibility. For example, near a magnetically ordered 
phase, such a susceptibility would reflect the enhance- 
ment of the magnetic fluctuations for k ~ Q, the an- 
tiferromagnetic ordering vector, and still more impor- 
tantly, the retardation effects implicit in the emergence 
of a new energy scale, "J" , associated with magnetic fluc- 
tuations. From this perspective, the various approaches 
to a "spin-fluctuation exchange" mechanism of supercon- 
ductivity appear as natural extrapolations of the present 
results to a more strongly coupled regime. However, it 
is important to stress that such procedures represent un- 
controlled approximations. There may well be competing 
channels or possibly a failure of the basic framework. One 
would like to have numerical calculations to test the util- 
ity of such phenomenological approximations. There are 
several specific studies we wish to discuss explicitly, as 
they have produced results which are particularly close 
to those obtained here: 

The two step RG approach of Schulai£ has already 
been mentioned. In the case considered there, the role 
of flo was played by a physical energy scale, defined as 
the energy scale below which the singular structures due 
to the proximate van Hove singularity could be ignored, 
and above which the system might as well be precisely 
tuned to the van Hove point. In particular, in that case, 
Oo (which they call t x )r^ takes on a U independent value 
which depends, instead, on t, t', and the chemical poten- 
tial, \x. Strictly speaking, unless the parameters are fine- 
tuned to be parametrically close to the van-Hove point, 
this physics does not survive the small U limit, as we have 
shown; our method includes, already, all effects of arbi- 
trary band-structure both at and away from the Fermi 



surface. 



VIII. DISCUSSION AND FURTHER 
DIRECTIONS 

The results we have obtained are asymptotically exact 
in the limit U/t — > 0, so long as the conventional RG 
treatment of the Fermi liquid is valid. In this limit, T c 
tends rapidly to zero, so the present results cannot be 
directly associated with a mechanism of "high temper- 
ature superconductivity." Moreover, it is clear that in 
most materials of interest, the interactions are moderate 
to strong, so the results cannot be said to have any direct 
relevance to these materials. However, some aspects of 
the results seem worth emphasizing which, when extrapo- 
lated to stronger coupling (where controlled calculations 
are not possible), may give insight into mechanisms of 
high temperature superconductivity in real materials. 

Most importantly, we have shown that repulsive in- 
teractions combined with lattice induced band-structure 
effects generically do result in a superconducting ground- 
state with non-trivial transformation properties with re- 
spect to the point group symmetries of the crystal. 

1) Where band-structure effects are weak, i.e. where 
the Fermi surface is nearly circular or spherical, the dom- 
inant superconducting instability is, generically, a two or 
three-fold degenerate spin-triplet p-wave. The driving 
force for this is more or less the same as originally envis- 
aged by Kohn and Luttinger. However, here if we com- 
pare T c with fixed small U/t at different band-fillings, 
the values of T c found in this regime arc small compared 
to the (still small) values obtained where the Fermi sur- 
face is more structured. Since in all the cases we studied, 
the p-wave state is a 2 or 3 fold degenerate representa- 
tion of the point group, and the d vector is also arbi- 
trary, there are many different possible forms of ordered 
state possible in principle. However, general consider- 
ations suggest that the generic ground-state will either 
break time-reversal symmetry, forming a "topological" 
p+ip superconductor (which corresponds to the A-phasc 
of Hclium-3 and presumably is the state that is observed 
in S^RuO^i), or the system will spontaneously break 
the relative "spin-orbit" symmetry as is the case in the 
B-phase of Helium-3. The perturbation which lifts the 
degeneracy between these two possibilities is the spin- 
orbit coupling of the normal state. Thus, an interesting 
extension of our work would be to incorporate the effects 
of spin-orbit coupling into the present analysis. 

2) On the square or tetragonal lattice with n near 1, 
higher values of T c are obtained in a spin singlet d-wave 
channel. This is encouragingly similar to what is found 
in the cuprate high temperature superconductors. 

3) On the triangular, hexagonal, cubic, FCC, and BCC 
lattices, there is a range of electron concentrations for 
which the dominant superconducting state is a two or 
three-fold degenerate singlet d-wave state. General ar- 
guments suggest that, in this case, the superconducting 
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state will be a time-reversal symmetry breaking d+id su- 
perconducting state. While we do not know of a currently 
well characterized material in which such a superconduct- 
ing state has been observed, the state is analogous to 
the "anyon superconducting state" originally proposed 
by Laughlini&i2. 

4) On both the triangular and hexagonal lattices, when 
there is more than one electron or hole pocket, the super- 
conducting ground-state is a non-degenerate spin triplet 
state, in which there is a full gap everywhere on the 
Fermi surface, but the gap changes sign in going from 
one pocket to the other. This situation gives the high- 
est transition temperatures (for fixed small U /t) that we 
have found. Moreover, the sign alternation on different 
pieces of the Fermi surface (although not the triplet char- 
acter) is reminiscent of the proposed gap function sym- 
metry in the Fe-pnictides. 

Conversely, it is interesting to study what aspects of 
the physics of real materials are inconsistent with the be- 
havior of the Hubbard model at small U/t. Taking the 
cupratc high temperature superconductors as an exam- 
ple, there are several features worth mentioning: 5) In 
the cuprates, T c is not very small in comparison to the 
microscopic energy scales. Indeed, it reaches values that 
are, within a factor of two or three, equal to the strong- 
coupling dimensional analysis estimate T c ~ |1 — n| J. 
In weak coupling, by contrast, T c is an emergent energy 
which is exponentially smaller than any of the micro- 
scopic energies in the problem. (However, T c does ex- 
trapolate to a value of order t in the limit pU ~ 1, as- 
suming that the extrapolation remains valid even where 
the justification for the result fails.) 6) In the cuprates, 
there is a clear breakdown of Fermi liquid theory at 
temperatures above T c , except possibly in the case of 
the most overdoped materials. Moreover, there is good 
evidence^— of substantial superconducting phase fluc- 
tuations persisting well above T c . In weak coupling, 
there is an emergent exponentially small crossover scale, 
flpT, at which perturbation theory breaks down, and 
where a corresponding breakdown of Fermi liquid the- 
ory is possible. However, because T c is exponentially 
smaller than Q pt , Fermi liquid behavior applies in a wide 
range of temperatures above T c - flpT 3> T > T c . BCS 
mean-field theory should, likewise, provide an extremely 
accurate description of the superconducting transition. 
In particular, the characteristic energy scale for phase 
fluctuations^ Tg = (n s /2m*)£<f ~ 2 - E F (£ /a) D - 2 , is 
exponentially larger than T c , precluding any substantial 
role for phase fluctuations. (Here n s /2m* and £o are, re- 
spectively, the zero temperature superfluid stiffness and 
coherence length, and a is a lattice constant.) 

3) Direct measurements of the gap function in the 
cuprate a 56 ' 57 suggest a gap function with the simple form, 
= [cos(k x ) — cos(ky)], at least in optimally doped 
materials. In real-space, this form implies the pair-field 
extends only to nearest-neighbor sites. In contrast, the 
weak-coupling gap function is more structured in k-spacc, 
as shown in Fig. [3] This reflects the non-local character 



of the induced pairing interactions and is, we believe, a 
generic feature of weak coupling. 

4) In the cuprates, superconductivity emerges from 
a doped antifcrromagnetic insulator with well defined 
spin-wave modes similar to those expected for a spin- 
1/2 Hcisenberg antifcrromagnct with exchange coupling 
J < t (and possibly some higher order exchange cou- 
pling representing^ the fact that U ~ 8t is not all that 
large). Moreover, there is direct evidence from neutron 
scattering that spin- wave-like collective excitations of the 
system with the same characteristic energy scal o 59 i 60 J 
persist into the superconducting phase at energy scales 
larger than the gap, even when there is no corresponding 
broken symmetry. In weak coupling, an antifcrromag- 
netic insulating state occurs only when two parameters 
are exponentially fine tuned, \t' /t\ < S and \(n — 1)| < S, 
where S ~ exp[ —a^Jt/U ]. A metallic antiferromag- 
netic state (as well as various other more exotic ordered 
phases) could occur with the fine tuning of only a sin- 
gle parameter, \(n — n v h)\ < $, in the vicinity of a van- 
Hove singularity. However, outside of this exponentially 
narrow range of concentrations, there is no identifiable 
local antiferromagnetic order present at weak coupling. 
Indeed, without some form of exponential fine-tuning, 
"competing orders" is not a concept that occurs in the 
weak coupling limit. 

There are many interesting directions in which the 
present work can be generalized. 

1) Most straightforwardly, it can apply to the Hubbard 
model on lattice systems with other band-structures than 
those considered here. 

2) It is clearly also possible to extend the same sort 
of analysis to situations in which there are more compli- 
cated interactions than the Hubbard U, such as multi- 
band models, where there are both intra-band and in- 
terband interactions, and even single band models with 
longer range interactions, such as a nearest-neighbor re- 
pulsion, V. Technically, these interactions complicate 
the analysis in the sense that a large number of other 
diagrams, involving the interaction between electrons of 
like spin, enter the pertubative analysis. Moreover, there 
are a variety of ways that the asymptotic analysis can be 
carried out, which will clearly lead to different physics. 
For instance, with interactions U and V, the analysis 
can be carried out as !7/t with U/V held fixed, or 
with U/V ~ (U/t) v where y is an appropriate positive 
exponent. 

3) Perhaps the most interesting extension involves the 
fine-tuning discussed above, which may permit the study 
of various strong-coupling phenomena in the weak cou- 
pling limit. In particular, the physics of competing orders 
can be explored by performing an asymptotic analysis in 
the limit U/t — > and n — > n v h in such a way that 
\n — n v h\ exp[ oty/l/U ] — > constant. We expect that 
this will permit us to reproduce much of the interesting 
physics obtained in various FRG calculations, but in a 
more controlled fashion. In order to study the interplay 
between superconductivity and a Mott insulating (an- 
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tifcrromagnctic) phase, an even more complex analysis 
must be performed, taking the limit U/t — > 0, n — > 1, 
and t'/t — > in such a way that |n — 1| exp[ a^/l/U ] — > 
constant and \t' /t\ exp[ a^l/U ] — > another constant. 
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Had we been interested in the negative U Hubbard model, 
r s (L a ) and T s (Li,) are the most important quantities and 
give rise to the BCS instability in the trivial s-wave chan- 
nel. Note that r(if,) eliminates the dependence on the 
initial cutoff fio as discussed in section VI. 

For unconventional superconductivity that derives 
from repulsive interactions, we will need to consider 
higher order diagrams. The next set of diagrams which 
are most important correspond to the set of diagrams in 
(L c ) of Fig. [ini Of these the second-order vertex shown 
in (2a) corresponds to 



r s (2a) = ^G (p)G (k + q+p) 



X (k + q) + O(n ) 



(A.5) 



Adding the two second-order terms, we find the expres- 
sion in Eq. [32] 

The most divergent third order term is given by dia- 
gram (Lrf) in Fig. [TH] which is simply 



Appendix: Perturbation theory in the Cooper 
channel 

The perturbative corrections to the 2-particle interac- 
tion vertices in the particle-particle channel are shown in 
Figs. H5l and [T6l In these diagrams two particles with mo- 
mentum k and — k are scattered to states with momenta q 
and — q. In the singlet channel, the incoming and outgo- 
ing pair of particles have the opposite spin whereas in the 
triplet channel, we require them to have the same spin. 
We first consider the diagrams in the singlet channel. 

The first order contribution to the interaction give in 
diagram (L a ) of Fig. \W\ is simply 



T s (L a ) 



1 



(Al) 



The first order term is independent of momentum and af- 
fects only the trivial s-wave pairing states. Next consider 
the higher-order terms. Before proceeding, we introduce 
as a notational convenience the following definition 



£ 



dp d d p 



n (27r) 



d+l 



(A.2) 



where po is a Matsubara frequency. The second-order 
ladder diagram shown in (Lf>) of Fig. [TBI is 



r s (L b ) = J2Go(p)g (-p) 



d d p 



2ep 



n (27r)' 

= P \og[A/n } + o(n ) (A.3) 

where Go(p) is the non- interacting Matsubara Green 
function 



G (p) 



1 



(A.4) 



T s {L d ) 



Y,Go(p)G (-p) 



= p 2 log 2 [ J 4/O ] (A.6) 



The next most important contribution comes from the 
diagrams in (L e ) and (Lj) with the thick solid line treated 
to second order (i.e. using diagram (2a)). These produce 

Tf\L e ) = J2 G ^+P + P r )Go(p r )G (p)Go(~p) 



X (q+p)[G Q (p)G Q (~p)]+O({l ) 
+ST s (L e ) (A.7) 



with 
ST s (L e ) 



PlWo) -x(q+P,fy] G {p)G {-p) 

(A.8) 

and the integrals in 6T s (L e ) are carried out with ilo = 
0, since this quantity is non-singular, and xiPlWo) is 
the frequency-dependent susceptibility. In practice, this 
quantity will be evaluated numerically. Upon summing 
over the Matsubara frequencies and making use of the 
identity 



d d p 
WT 



dp vf 
Sf v F {p) 



dep(e) 



(A.9) 



|e|>O 



where e is the energy relative to the Fermi surface, it 
follows that 

ri 3 )(L e ) = 7 (3) (<?)piog [A/a ] + sr s (L e ) (A.io) 

The remaining third order terms are non-singular and arc 
obtained from diagram (L c ) again, but this time, keeping 
only 3rd order contributions to the thick solid line, which 
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are obtained from diagrams (3a)-(3g). While diagram 
(3a) is simply 

T s (3a) = Y, G o(p) G o(k + q+p)G (p')G Q (k + q+') 
pp' 

= X 2 (k + q)+O(n ), (A.ll) 
diagrams (3b, 3c) which are 

T s (36) = F s (3c) = J2 G a (p)G a (q+ P '-p)Go(p')G Q (k+q+p') 
pp' 

(A.12) 

are not expressible in simple closed form due to the in- 
ternal frequency integration (but it is possible to ob- 
tain these quantities to arbitrary precision in a numerical 
computation). Diagram (3d) contributes 

T s (3d) =x 2 (k-Q) (A.13) 

which is also non-singular, and lastly, the diagram of the 
form given in (3e) contributes 

T s (3e) = - Go(p)G (p + k - q)Go( P ')G (p' -k-p) 
pp' 

(A.14) 

Next, we consider the fourth-order terms. Again, the 
most singular contribution comes from the ladder shown 
in (L g ) of Fig. M 

T(L g ) = p 3 log 3 [A/too] + O(fio) (A.15) 

The next most important contributions to this order 
comes from the diagrams in (L^,Li) with the thick solid 
line treated again only to second order, namely from di- 
agram (2a). Diagrams (L^) and (Lj) have the following 
contributions to leading logarithmic order: 

r s (L /i )=p 2 log 2 L4/Q ] 7 ( 3 >(fc) 
r s (L i ) = p 2 \og 2 [A/Sl ] 1 W(q) (A.16) 

where again, we have neglected the frequency dependence 
of the susceptibility (which does not produce a singular 
contribution), and have made use of the identity in Eq. 
IA.9I Note that F S (L/,) depends only on the incoming 
momenta, whereas T s (Li) depends only on the outgoing 
momenta. In this way, it is easy to see that diagram (Lj) 
has no momentum dependence whatsoever: 

L S (L,) = p 2 log 2 L4/f! ] 7l 4) (A. 17) 



The next leading terms are obtained again from the di- 
agrams (L e , Ly), except now, we require that the thick 
black lines in these diagrams take the non-singular third- 
order vertex which is given by Tg described above: 

if >(v) = P \og[A/n,]^\q) + o(n Q ) (A.i8) 

and a similar expression holds for T^iLf). The final 
term which is singular at fourth order is given by the di- 
agram (Lfc), where the thick solid line is approximated by 
the second-order contribution in (2a). This contributes 

T s (L k ) = Y,x(k + P)G a (p)G (-p)x(p + q) 
p 

= U 4 plog [A/n ] 7 i 4) (fc, 0) + O(O )(A.19) 

The remainder of terms that occur to fourth order do not 
make a singular contribution to the vertex. We shall not 
discuss them here. 

Now, we consider the triplet channel which has con- 
siderably fewer diagrams. The lowest order contribution 
comes from the second-order diagram in (2b): 

T t (2b) = -x(k-q) + O(il ) (A.20) 

The third order correction to the vertex comes from dia- 
gram (3g) which is 

r t (3<?) = -Y / Go(p)Go(k-q + p)Go(p')G (k+p + p / ) 
pp' 

= - £ Go(p)G (k - q + p) x (k + P ) + 0(^)21) 

p 

Lastly, the fourth-order term comes from diagram (L„), 
treating the thick solid line to 0(U 2 ), ie. using diagram 
(2b). This produces 

T t (L n ) =p\og[A/n } 1 { t 4} (k,q) + O(n ) (A.22) 

This contribution, although higher order, has an impor- 
tant conceptual importance: it has the correct logarith- 
mic dependence on fJ to remove the dependence of the 
final expression for T c on the initial cutoff. 
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